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A measure of the impact of future dark energy experiments based on discriminating 

power among quintessence models. 
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We evaluate the ability of future data sets to discriminate among different quintessence dark 
energy models. This approach gives an alternative (and complementary) measure for assessing 
the impact of future experiments, as compared with the large body of literature that compares 
experiments in abstract parameter spaces (such as the well-known wo — w a parameters) and more 
recent work that evaluates the constraining power of experiments on individual parameter spaces of 
specific quintessence models. We use the Dark Energy Task Force (DETF) models of future data 
sets and compare the discriminative power of experiments designated by the DETF as Stages 2, 
3, and 4 (denoting increasing capabilities). Our work reveals a minimal increase in discriminating 
power when comparing Stage 3 to Stage 2, but a very striking increase in discriminating power when 
going to Stage 4 (including the possibility of completely eliminating some quintessence models). We 
also see evidence that even modest improvements over DETF Stage 4 (which many believe are 
realistic) could result in even more dramatic discriminating power among quintessence dark energy 
models. We develop and demonstrate the technique of using the independently measured modes of 
the equation of state (derived from principle component analysis) as a common parameter space in 
which to compare the different quintessence models, and we argue that this technique is a powerful 
one. We use the PNGB, Exponential, Albrecht-Skordis, and Inverse Tracker (or Inverse Power 
Law) quintessence models for this work. One of our main results is that the goal of discriminating 
among these models sets a concrete measure on the capabilities of future dark energy experiments. 
Experiments have to be somewhat better than DETF Stage 4 simulated experiments to fully meet 
this goal. 
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I. INTRODUCTION 

Over the last decade or so there has been mounting 
evidence that the energy of the universe has a large ac- 
celerating component, dubbed "dark energy" [l|. As the 
evidence has become more convincing, there has been 
growing enthusiasm for launching a major program to 
collect additional data that will help us better under- 
stand the nature of the dark energy 0] , and indeed con- 
siderable progress is being made on this front. 

To plan a strong program of dark energy studies one 
needs to assess the relative impact of different possible 
experiments. This has most often been done by describ- 
ing the dark energy in some abstract parameter space 
and calculating how much a given data set could con- 
strain those abstract parameters. For example, the Dark 
Energy Task Force (DETF) [3J, building on earlier work 
0, Q , used a standard two parameter model of the dark 
energy equation of state w as a function of cosmic scale 
factor a given by w(a) — wo+w a (l — a) to form a figure of 
merit based on constraining power in the wq — w a param- 
eter space. For the most part, other authors have used 
other abstract dark energy parameterizations [3, [g, Qj], 
but more recently we have extensively explored the im- 
pact of future experiments using the constraints pro- 
duced on the actual parameters of scalar field dark energy 
models [8|, [9|, LLCJ, LU| • That work gives another window on 
the power of future experiments, which we have argued is 
largely consistent with the DETF results in the wq — w a 
parameter space. 



Our recent work [8|, [jj, USi, LUj shows the constrain- 
ing power of future experiments on specific dark energy 
models. However, because the natural parameter space of 
each quintessence model is very different from the oth- 
ers, we were not able to use our techniques to directly 
evaluate the ability of experiments to favor one dark en- 
ergy model strongly over another. The one exception to 
this is "cosmological constant" dark energy (which has 
it; (a) = —1). Each quintessence model we considered had 
a part of parameter space where the quintessence closely 
mimicked a cosmological constant, and we used that fact 
systematically in [8|, [tj QiJ, to consider discriminat- 
ing power in the quintessence vs. cosmological constant 
domain. 

This paper builds on that earlier work to make a more 
comprehensive analysis of the ability of future data to 
discriminate among different quintessence models. The 
key new ingredient is the use of a specially chosen param- 
eter space to represent the different quintessence mod- 
els in a common and comparable form. To this end 
we use the "independently measured modes" of w(a), 
which have long been appreciated for a variety of rea- 
sons @, 0, [H Gl Q M, US G3- The modes, which 
are different for each experiment, represent the modes 
or "moments" of w(a) of which uncorrelated measure- 
ments are made by that particular experiment. This fea- 
ture allows us to identify the modes which are best mea- 
sured and analyze them in a straightforward way (due 
to the lack of correlations) . These modes comprise a ba- 
sis which spans the space of possible functions w(a). A 



given quintessence model with specific fixed parameters 
will give a specific function w(a) which can then be ex- 
panded in the modes, and the expansion coefficients form 
the parameter space in which we work. 

We use the PNGB [|, [3, Exponential 0, El, 
Albrecht-Skordis Jal2IJ|' anc ^ Inverse Tracker (or Inverse 
Power Law) [ll|, |21| quintessence models for this work. 
This is a diverse set of models each of which holds its 
own special interest among researchers (see our discus- 
sion in the introductions of J8L |fj, LLQl LLJJ for a brief re- 
view of the motivations and [22| to place these models in 
a more general context). One of our results is that these 
four quintessence models actually occupy very different 
regions of the mode parameter space. This tells us that 
the aspects of w(a) that are well measured by realistic 
experiments have the potential to be extremely useful in 
discriminating among quintessence models. How much 
this potential is realized is of course related to the res- 
olution achieved by a given experiment in its mode pa- 
rameter space, and that issue is the subject of much of 
our analysis. 

One of our key results is that the goal of discriminating 
among these four models sets a concrete measure on the 
capabilities of future dark energy experiments. Experi- 
ments have to be somewhat better than DETF Stage 4 
simulated experiments to fully meet this goal. 

Sections ILT1 IIIIl and IIVI describe our methods and Sec- 
tion|V]presents our detailed results, while Section [Vl] pro- 
vides a discussion of the relevant statistical issues. Our 
conclusion, Section lYlIl summarizes our results. 



II. CONNECTION TO OUR EARLIER WORK 



Our work builds very directly on our recent papers 
studying specific quintessence models @, IjJ [l(J, [Tl[ . We 
refer the reader to those papers to learn more about 
our methods. (An appendix giving the technical details 
about our methods that are common to all these papers 
can be found in [9(.) One product of this earlier work is a 
set of Monte Carlo Markov Chains representing the dis- 
tribution of models that are consistent with specific Stage 
2 simulated data that is chosen to be consistent with 
a cosmological constant cosmology. Specifically, these 
chains represent the distribution of possible scalar field 
parameters that are consistent with a specific simulated 
Stage 2 data set. Each quintessence model has its own 
scalar field parameters and its own chain representing 
the distribution in that space. Also, in each case we base 
the Stage 2 data around a specific "fiducial" set of scalar 
field parameters that are consistent with a cosmological 
constant at the 1 — 2tr level. Using these different fiducial 
models accounts in a rough way for uncertainties in the 
outcomes of the Stage 2 experiments. 

We use these chains for the work reported here by de- 
termining the equation of state function w(a) for each 
point on the chain and projecting that function into 
the eigenmode-based space corresponding to a particu- 



lar simulated data set from Stage 3 or Stage 4 (as dis- 
cussed in detail below). When data from different scalar 
field models is analyzed using the same eigenmodes, those 
modes provide a common parameter space in which scalar 
field models with different "fundamental" scalar field pa- 
rameters can be compared on a common footing. In this 
usage the full chain represents how scalar field models 
that are consistent with each other at Stage 2 would be 
distributed in the eigenmode-based space based on data 
from Stage 3 or 4. Any discriminating power among re- 
gions of the Stage 2 chains enabled by the higher stages 
represents progress over Stage 2 data, on which the orig- 
inal chains are based. 



III. GENERATING EIGENMODES 

For this work, eigenmodes were generated using 
MCMC calculations. We do this by breaking up the 
equation of state into nine bins linear in scale factor from 
a = .2 to a = 1 and using the value of w + 1 in each of 
these bins as parameters. We then run Markov chains 
with these "bin" parameters in addition to cosmologi- 
cal parameters in order to calculate a covariance matrix, 
from which we extract the nine-by-nine covariance ma- 
trix for said "bin" parameters. The eigenvectors of this 
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FIG. 1: The first three best measured eigenmodes of w(a). 
The Stage 3 (dot-dash), Stage 4 ground (dash), and Stage 
4 space (solid) modes share a common general form at each 
mode level. Specific differences can be related to various dif- 
ferences among the experiments including how deeply a given 
experiment probes in redshift. (Technically, these modes 
should each be represented by nine discrete bin values. The 
connecting lines guide the eye, and reflect a likely "continuum 
limit" as discussed in Q.) 



matrix give us our eigenmodes, while the eigenvalues are 
the squares of the uncorrected error in each mode. While 
these methods are slightly different from the Fisher ma- 
trix techniques of |7|, the results are consistent, and our 
choice of binning is driven by the analysis in [7]]. As 
with our previous work [S, M, LLtJ, LLU] , the simulated data 
sets are constructed in a manner equivalent to the DETF 
simulated data. We do not include cluster data (due to 
the technical difficulty of including it discussed in 0). 
There are a number of possible considerations beyond the 
DETF work (such as more carefully considering the im- 
pact of including cross correlations among different pho- 
tometric data types [23|, (24|) that many expect will lead 
to significant improvements over the DETF projections. 
For this paper we use the original DETF data models for 
ease of comparison, except briefly in section [V] where we 
consider a simpleminded extension. 

Plots of the first three eigenmodes are given in Figure 
[1] (ranked by how well each mode is measured). All nine 
eigenmodes together form an orthonormal basis, which is 
different for each data set. The modes pick up additional 
oscillations as one goes from best measured (mode f ) to 
less well measured modes. 



IV. PROJECTION ONTO EIGENMODES 

We use the above eigenmodes to analyze four 
quintessence models of dark energy: the Pseudo-Nambu- 
Goldstone Boson (PNGB), Exponential, Albrecht- 
Skordis (AS), and Inverse Tracker (IT) models. Sample 
equation of state behavior of these models is illustrated in 
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FIG. 2: A characteristic w(a) function for each of the four 
scalar field models considered in this paper. 



FIG. 3: Stage 3 photometric: A plot of the projections of 
the PNGB, Exponential, AS, and Inverse Tracker models onto 
the first four Stage 3 photometric eigenmodes. The displayed 
points are from MCMC chains for each model, and the scale 
of each axis is given by at, the uncertainty in measurements 
of that mode. 



Figure^] .26|. In each case, we will use points pulled from 
MCMC chains as representations of each model's param- 
eter space. The chains were run on DETF Stage 2 type 
data generated using a fiducial point in each model, as in 
[g, \Q, [l(| ; this gives us a fairly wide spread of parameter 
space for each model and represents in a rough way the 
uncertainties in the outcomes of the Stage 2 experiments. 
We use a simple algorithm to average the equation of 
state w(a) in the nine scale factor bins, Wj. These Wj 
can be mapped into mode projections rrii by matrix mul- 
tiplication: 



Y,E lJ (w J + l) = 



(1) 



where Ey is the j term of the i eigenmode. We use 
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FIG. 4: Stage 4 ground: A plot of the projections of the 
PNGB, Exponential, AS, and Inverse Tracker models onto the 
first four Stage 4 ground eigenmodes. While the first mode 
is measured better by this data than by the Stage 3 data, it 
is the improvement in the measurement of the second mode 
makes the models distinguishable over a wide range of their 
parameters. 



FIG. 5: Stage 4 space: A plot of the projections of the 
PNGB, Exponential, AS, and Inverse Tracker models onto 
the first four Stage 4 space eigenmodes. As with the Stage 
4 ground plot, the most significant improvement, in terms of 
model distinction, going from the Stage 3 data to the Stage 
4 space data is the measurement of the second mode. 



Wj + 1 to center the eigenmodes around w(a) = — 1, so 
that rrii = for all i is a cosmological constant cosmol- 
ogy. The actual value of rrii for a given w(a) depends 
on the number of bins used, so it is more convenient to 
look at — , where U{ is the square root of the i th covari- 
ance matrix eigenvalue. This expresses the power in each 
eigenmode in units of its uncorrelated error, and should 
be relatively stable as we change the number of modes 
by, for example, refining the bin size. 

When we consider how to display the range of 
quintessence models in the eigenmode space, graphing 
the power in the first three modes in a rotatable three 
dimensional plot can be fascinating. However, as this 
does not lend itself to the static two dimensions of a pa- 



per, it is more enlightening to view the modes two at 
a time (an interested reader may contact the authors of 
this paper to view the three-dimensional versions). 

Figure [3] shows Stage 2 distributions of our four ex- 
ample quintessence models represented in the first four 
modes of combined Stage 3 data (using photometric data 
and systematics designated "optimistic" by the DETF, 
just as was done in [8|, 0, [Hi EH)- Because the axes of 
these plots are scaled in units of the error of each mode, 
one can get an intuitive idea of the Stage 3 resolving 
power by noting that Stage 3 data should roughly resolve 
areas of unit size in these plots. 

It is interesting to note how the models examined here 
occupy very distinct portions of the "mode space" ex- 
cept for where they meet at the origin (w(a) = —1). 



We should mention here that these figures in many cases 
display only a subset of the total space accessible by the 
models, because the parameter spaces of the models were 
restricted in the MCMC chains. This is particularly true 
of the AS model, for which we expect that the kinks and 
wiggles in the displayed distribution will characterize a 
broader distribution of such features. This may also be 
an issue for the PNGB model, which was restricted in 
the MCMC analysis to the concave down portion of its 
potential. The Exponential model, however, has a con- 
cave up potential and gives us some idea of where the 
concave up portion of the PNGB model may lie in the 
mode space. As discussed in [3, M, UM, HH, m °st of these 
restrictions were required to eliminate parameters that 
would be completely unconstrained by even the best fu- 
ture data. 

On a similar note, there is a small fraction of the points 
(< 1%) in the Inverse Tracker model that never display 
tracking behavior, but instead display the thawing 25] 
behavior of the PNGB and Exponential models; this can 
be seen in the plots as a handful of outlying points above 
the main concentration of Inverse Tracker model points. 

It is of interest that, as we look at parameterizations 
of increasing distinction from the cosmological constant, 
we see a consistent increase in the amplitude of their 
first mode projections. As such, we might look at the 
first mode projection as a signal of the presence of one of 
these scalar field models, but only in extreme cases will 
it help us distinguish between them (noting that a very 
large value would actually rule out the AS model). It 
is clear that it is the higher modes that will distinguish 
between models, though the plot shows that Stage 3 data 
will not do this very well. 

The authors of [y| also made plots using similar modes 
to ours, but their approach differs. In [6( the potentials 
are drawn from an abstract continuous space that covers 
large ranges of possibilities, whereas our potentials are 
drawn from specific scalar field models which each only 
have certain classes of behaviors even as the parameters 
are varied fully. These differences have allowed us to 
discover a much more striking structure than can be seen 
in@. 

In Figure[5]and[5]we show the mode projection plots for 
Stage 4 space and ground data sets, respectively (using 
the same DETF "optimistic" combinations used in [g, [9|, 
flG . Ill|). We can see from these plots that the overall 
shape of the distribution of projected models remains 
very similar in the first two modes. The Stage 4 data sets 
do a better job of measuring the first mode than Stage 
3, a much better job of measuring the second mode, and 
even begin to give resolution of the third mode that is 
of marginal value in resolving these models. As a result, 
for a scenario in which the Stage 3 projects detected the 
presence of quintessence at the level of a few sigma, the 
Stage 4 data should be able to discriminate between this 
sample of models to at least that level. 
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FIG. 6: A plot of the mean (with bars at one standard devia- 
tion) and maximum projections of the examined models into 
the Stage 4 space mode-space, in terms of the mode uncer- 
tainties. It can be seen that for nearly all of the model points, 
the projections onto the modes higher than three are negligi- 
ble. Though the some models have outlying points that have 
some power in the higher modes, the bulk of each model does 
not. 



RESOLVING THE QUINTESSENCE 
MODELS 



With the models projected onto the eigenmode space, 
we then have a common space in which to compare the 
different scalar field models and evaluate how well com- 
ing experiments will distinguish among them. We can 
use mode uncertainties as a metric to find how measur- 
ably different each quintessence model will be. Informed 
by the plots of mode projections, we can reasonably ex- 
pect that the first mode will dominate this measure, but 
the second and third modes will play a part in setting the 
minimum separation. As for the higher modes, we can 
expect to see these swamped by their uncertainty; a com- 
bination of the poor measurement and the lack of power 
in these modes by our quintessence models makes them 
largely irrelevant. Figure [6] gives information about the 
distribution of mode amplitudes for our chains. From 
Fig. [6] one can see that the measurement error in the 
higher modes will be much larger than the distribution 
of physically interesting values. This gives us, in effect, a 
prior on the value of these modes that is much stronger 
than the data constraints from even Stage 4 experiments. 
Therefore, this calculation will actually be done using 
only the first three eigenmodes. 

Assuming Stage 2 data as discussed above, higher- 
stage data might come from a universe based on any 
of the quintessence models represented by points in our 
mode space plots. Our next step is to scan all possi- 
ble data outcomes and evaluate the potential discrimi- 



2 
1 


to - 1 

CM _ 2 

CD 

o 3 

E _ 4 

-5- 
-6- 

-7 



• PNGB 

+ Exp. 

X Inv. Trac. 

O AS 



O 

o 
o 



o 



6 8 10 12 14 

mode 1/c 



~5 

CD 
O 



X 



PNGB 
Exp. 
Inv. Trac. 



O AS 



X 



O x 



X 

O 



O 



O 



0.5 1 1.5 2 2.5 3 

mode 3/cr 

FIG. 7: The four test points for each model (16 test points 
in all) represent possible experimental outcomes. The test 
points are plotted here in the eigenmodes for Stage 4 space. 
One can compare this figure with Fig. [5] to see the degree to 
which the full range of behavior of each model is represented 
by the test points. 



nating power. We represent the Stage 2 distribution of 
each quintessence model, by choosing four "test points" 
spread evenly across the mode space. Four test points 
each for four quintessence models gives a total of six- 
teen "data outcomes" which are meant to represent (in a 
rough way) the full range of possibilities. The distribu- 
tion of our test points in mode space is shown in Fig. [7] 
These "test points" were chosen by their first mode pro- 
jections, which places them in sequence along the very 
nearly linear regions the models cover [27j |. For each of 
the sixteen test points, we analyze the ability of data 
centered on the test point to exclude other points on the 
chains. Note that in Fig. [7] the test points are not rep- 
resented with noise, which is expected in any data set 
(at a level given by the a' t s) and which is reflected in our 
quantitative analysis. 

The curves in Figs. El [9j and [IJJ show the fraction 
of model points with x 2 less than a specific value, given 



TABLE I: Stage 3 photometric. This is a table of the min- 
imum x 2 (ignoring the smallest 1% of x 2 values). The 99% 
confidence level is a x 2 of 11.36. Notably, the only compar- 
isons that rise to that level are between the test points with 
large first mode projections and the AS model, telling us that 
this is mostly a first mode measurement. As discussed in the 



text, this table gives the values of x 
[8] sharply approach the x-axis. 



where the curves in Fig. 
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on the a;-axis[28|. To read these plots, start with the 
labels on the left hand side. These read PNGB (Pseudo- 
Nambu-Goldstone Boson), EXP (Exponential), IT (In- 
verse Tracker), and AS (Albrecht-Skordis), and label 
which model the "test points" used in that row of plots 
were pulled from. Then look to the top of the figure, 
where the same model labels mark the columns; these 
label which model the test points are being compared to 
in each column. In each plot, four functions are graphed, 
X 2 vs. the fraction of the compared model's points that 
have that x 2 or lower relative to the test point in ques- 
tion. For example, looking at the PNGB vs. PNGB plot, 
the curves represent how far the rest of the PNGB model 
is from each PNGB test point. The point where each 
function touches the x-axis (the minimum x 2 ) is in this 
case a loose measure of how densely the MCMC chain 
populates the mode space at that test point. However, 
if one were to look at the PNGB row and the AS col- 
umn, the graph there represents how far the PNGB test 
points are from all of the points in the AS chain, and 
the minimum x 2 of each function represents by what x 2 
each of the PNGB test points would rule out the entire 
AS model. 

The sharp left-hand edges of curves in Figs. [HI 03 and 
1101 imply very strong discriminating power at the level 
of x 2 given by the point on the x-axis where the edge is 
located. We have organized information about this im- 
portant feature in the following tables: Tables [H HH and 
IIIII give numerical results for the minimum x 2 one could 



TABLE II: Stage 4 Ground. This is a table of the minimum 
X 2 (ignoring the smallest 1% of \ 2 values). Again, the low 
first mode projections for the AS model are responsible for the 
highest x 2 values, but we can also see significant separation 
for test points that have first mode projections in the range 
of the model they are being compared to. As discussed in the 



text, this table gives the values of \ 
[9] sharply approach the x-axis. 



where the curves in Fig. 
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TABLE III: Stage 4 space. This is a table of the minimum x 2 
(ignoring the smallest 1% of x 2 values). Again, the low first 
mode projections for the AS model are responsible for the 
highest x 2 values, but we can also see significant separation 
for test points that have first mode projections in the range 
of the model they are being compared to. As discussed in the 



text, this table gives the values of x 
[10] sharply approach the x-axis. 



where the curves in Fig. 
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estimate from Figures [8] [9] and [10] In general (due to 
the outlying points and other factors), one expects some 
low level tails even on the otherwise sharp rising edges. 
In order to make sure we are quantifying the true edge 
of the figure we ignore the closest 1% of points from the 
compared model. This is equivalent to finding the x 2 
where each plot crosses the .01 fraction mark. For com- 
pleteness, we have again included the Exponential model, 
though we do not expect any experiment to distinguish 
between it and the PNGB model. The numbers reported 
here reinforce the intuition we gain from the mode plots: 
only in extreme cases will Stage 3 distinguish between 
the models. The only examples that rise above the level 
of 99% confidence (which is a x 2 of 11.36 for the three pa- 
rameters used here) are due to the previously mentioned 
observation that the AS model has a much smaller range 
of amplitudes in the first mode than the other three mod- 
els. 

The Stage 4 ground and Stage 4 space x 2 values show 
significant improvement over Stage 3. Again, the largest 
X 2 values come from the large first mode separation be- 
tween the last test point for each model and the nearest 
AS model point. 

Looking at Tables HTl and lllll it appears that the Stage 4 
data lie on some kind of threshold: There are quite a few 
entries greater than 11.36 (indicating exclusion at 99% 
confidence) but plenty that are lower. To explore the na- 
ture of this threshold a bit more, Table IIVI presents x 2 
values for a hypothetical data set that would improve the 



TABLE IV: Stage 4 space, with experimental uncertainty re- 
duced by 2/3 in each mode. This is a table of the minimum 
X 2 (ignoring the smallest 1% of x 2 values). 
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(Ji by a factor of 1.5 over Stage 4 space. Due to a variety 
of considerations (including those discussed in [23|, l24j ) 
many believe such improvements (or even much better 
ones) over the DETF modeling of Stage 4 to be realistic 
for some experiments. In Table IIVI there are indeed a 
great many more entries greater than 11.36, further sup- 
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porting the notion of a threshold around Stage 4. As a 
reference point, the ctj of the Stage 4 data sets are a fac- 
tor of about 4 or 5 smaller than the Stage 3 o^ for most 
modes. To get a more complete understanding of which 
experiments might achieve (or beat) the level of model 
discrimination shown in Table HVl one would have to un- 
dertake detailed modeling of alternative experiments and 
data analysis schemes. Such a systematic analysis is not 
the subject of this paper, and we leave that for future 
work. 



VI. DISCUSSION OF STATISTICS 

The MCMC chains of the selected scalar field models 
in this analysis are used as a means to represent the phe- 
nomenology associated with each model, that is, to rep- 
resent the full spread of points in mode space the model 
can occupy. This is not to be confused with the set of 
priors that would be used in a Bayesian analysis. The 
question we seek to answer with Figs. 13 E3 and[lOl and 
Tables HI to II VI is "If the universe is a single realization of 
a particular model, how well will a given experiment rule 
out other models/phenomenologies?" The method we 
use to answer this is to take the probability of observing 
the observables generated by that realization (for a given 
experiment) at all other points in the eigenmode space, 
and then find the highest of those probabilities (small- 
est x 2 ) ou t of our sets of model-based phenomenologies. 
We can then speak to the power of the experiments to 
distinguish between models. 

Because we are only interested in the regions of closest 
approach between a model and a hypothetical realization, 
the exact distribution of the points in the chains used for 
this work is largely unimportant, so long as it does a 
reasonable job of exploring the possibilities of each. Our 
focus on the regions of closest approach is really a luxury 
which is enabled by the remarkable structure we have 
found in the eigenmode space. An outcome with more 
overlap (such as we do see between the Exponential and 
PNGB models) might draw one into a Bayesian analysis 
where priors on the different model distributions in mode 
space would play a critical role. One can loosely think 
of our approach as a very conservative one that allows 
for very unprincipled theorists who might cook up rea- 
sons to place delta function priors on the point of closest 
approach. We have shown that high quality Stage 4 ex- 
periments could rule out entire models, even under these 
adverse circumstances. 

There is an important assumption built into our anal- 
ysis. Our choice of four test points spread out along the 
mode curves for each model implies that we are giving 
"equal weight" to the actual Universe residing anywhere 
along the mode space curves. In the current confused 
state of dark energy theory, one can expect a wide variety 
of individual preferences on this point. For someone with 
a different preference, our results about the structure in 
mode space remain unchanged, and will still impact any 



alternative analysis performed to express different preju- 
dices about possible outcomes of future experiments. We 
feel we have built assumptions into our analysis which 
allow us to most directly address the question "if one of 
these experiments is carried out and gets a significant 
signal, what are the prospects for that signal completely 
ruling out a given scalar field model?". While we have 
focused our quantitative analysis on cosmologies based 
on one or the other of our scalar field models, one also 
can see (by inspecting the mode space figures) that it 
would be quite possible to get an experimental signal 
that rules out all the scalar field models considered here. 
That again is a consequence of the mode space structure 
demonstrated in this paper. 



VII. CONCLUSIONS 

We have considered the ability of future dark energy 
experiments to discriminate among different scalar field 
quintessence models of dark energy. To this end we have 
projected the equation of state functions w(a) for each 
model into the space of best measured eigenmodes of fu- 
ture experiments. We believe this approach is effective 
and convenient for investigating the ability of a given 
data set to discriminate among different quintessence 
models. Specifically, this approach offers a way around 
the fact that parameters of different quintessence models 
are typically not defined in the same spaces, which makes 
more direct comparisons of the models problematic. 

The four quintessence models considered here cre- 
ate a distinctive structure when projected into the 
mode spaces. The goal of discriminating among these 
quintessence models gives an alternative and complemen- 
tary measure of the impact of future experiments. In 
large part due to the structure in mode space, this mea- 
sure has some striking features that are different from 
other measures considered previously. 

We have shown that the DETF Stage 3 data will 
have very little utility in discriminating among the four 
quintessence models, although it will significantly probe 
the possibility of non-cosmological-constant-like behav- 
ior. DETF Stage 4 simulated data appears to lie right at 
an interesting threshold in that this data shows signifi- 
cant discriminating power among the quintessence mod- 
els we considered. We also showed that modest improve- 
ments over DETF Stage 4 (which many consider quite re- 
alistic for some Stage 4 experiments currently proposed) 
allows one to cross this threshold more completely, lead- 
ing to substantially greater discriminating power. 

It is important to note that at our current level of theo- 
retical understanding all quintessence models are suspect, 
and we are not advocating the use of the measures pre- 
sented here to the exclusion of other approaches. How- 
ever, as discussed in y§, M, [[(J, UM > we hare chosen an in- 
teresting sampling of reasonably motivated quintessence 
models. Since such quintessence models are part of the 
current theoretical discussion of dark energy, discrimi- 



nating power among these models should be part of how 
we evaluate the impact of dark energy experiments. 
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FIG. 8: Stage 3 Photometric. These plots show the distributions of x as each test point is compared with each comparison 
model. Specifically we plot the fraction of \ 2 values less than \ 2 given on the :r-axis. The rows of this figure correspond to 
which model the test points have been pulled from (as seen in Figure [7}, while the columns correspond to the model which is 
compared to those test points. In each plot, the test point closest to the origin of the mode space (and thus, a ACDM model) 
is denoted by a dotted line, the next closest as a dot-dashed line, the next as a dashed line, and the farthest as a solid line. 
The 99% confidence interval for three parameters is \ 2 = 11-36. The relatively sharp left-hand edges of these curves are an 
interesting feature (related to the gaps between models in mode space) which is discussed in the text and Table [I] 
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FIG. 9: Stage 4 Ground. These plots show the distributions of \ 2 as each test point is compared with each comparison 
model. Specifically we plot the fraction of x 2 values less than \ 2 given on the a;-axis. The rows of this figure correspond to 
which model the test points have been pulled from (as seen in Figure [7}, while the columns correspond to the model which is 
compared to those test points. In each plot, the test point closest to the origin of the mode space is denoted by a dotted line, 
the next closest as a dot-dashed line, the next as a dashed line, and the farthest as a solid line. The 99% confidence interval 
for three parameters is \ 2 = H-36. The relatively sharp left-hand edges of these curves are an interesting feature (related to 
the gaps between models in mode space) which is discussed in the text and Table ITT1 
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FIG. 10: Stage 4 Space. These plots show the distributions of % 2 as each test point is compared with each comparison model. 
Specifically we plot the fraction of x 2 values less than \ 2 given on the :r-axis. The rows of this figure correspond to which model 
the test points have been pulled from (as seen in Figure [7]), while the columns correspond to the model which is compared 
to those test points. In each plot, the test point closest to the origin of the mode space is denoted by a dotted line, the next 
closest as a dot-dashed line, the next as a dashed line, and the farthest as a solid line. The 99% confidence interval for three 
parameters is % 2 = 11-36. The relatively sharp left-hand edges of these curves are an interesting feature (related to the gaps 
between models in mode space) which is discussed in the text and Table Hill 



